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ABSTRACT 


Measurements of electron densities in the D region are made by the 
partial-reflection differential-absorption and differential -phase techniques. 
The differential-phase data are obtained by a hard-wired phase-measuring 
system. Electron- density profiles obtained by the two techniques on six 
occasions are plotted and compared. Electron-density profiles obtained at 
the same time on 30 occasions during the years 1975 through 1977 are averaged 
to form a single profile for each technique. The effect of .varying the 
assumed coll is ion- frequency profile on these averaged profiles is studied. 

It is found that adjusting the assumed collision frequency will bring the 
averaged electron-density profiles into close agreement, suggesting an error 
in the collision frequency. 

Time series of P-region electron- density data obtained by 3.4-minute 
intervals on six days during the summer of ' 1977 are examined for wave- like 
disturbances and tidal oscillations. Short- and long-term oscillations are 
observed, and the coherency spectra of these oscillations indicate that some 
of these fluctuations are possibly due to propagating internal gravity waves. 
The diurnal variations on the six days show an asymmetry about local noon, 
with a bite-out at noon. Plots of daily electron densities at 72, 76.5, and 
81 km for the period of October 1975 through April 1976 show increased vari- 
ability during the winter months, a clear indication of the occurrence of 
» 

the winter anomaly at Urbana. 
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1. INTRODUCTION 

1 . 1" D-Region Characteristics 

The D region is the lower part of the ionosphere, extending from 60 to 
90 km. The D region is the least understood, and most complex region of the 
ionosphere. The daytime upper D - region ionization is due mainly to photo- 
ionization of nitric oxide by solar Lyman-a radiation [Strobelj 1971] . Other 
sources of ionization include solar X-radiation of major D-region constituents 
(Og and N 2 ) , galactic cosmic radiation of constituents in the lower D region 
(60 to 70 km) , and extreme ultraviolet radiation of metastable molecular 
oxygen, 02(*A) [Whitten and Poppoff, 1971] . Also, precipitating energetic- 
electrons may be a .significant source of ionization at midlatitudes during 
geomagnetic storms. The relatively high neutral atmospheric pressure is 
responsible for the complexity of the ion chemistry in the D region. Ef- 
ficient three-body reactions may occur and large water cluster ions may form. 

In the lower D region the relatively high neutral density allows rapid three- 
body attachment of electrons to molecular oxygen, which reacts with other 
neutral gases to form other negative ions. 

Because the D region is formed mainly by solar radiation, it almost 
completely disappears at night. The average values of electron densities are 
higher in the summer months due to a lower solar zenith angle. Electron 
densities show many variations due to dynamical effects. Variations are 
observed on many time scales and are caused by many processes. In this work 
seasonal, day-to-day, diurnal, and shorter time scale variations will be' 
studied. 

✓ 

1.1.1 The winter anomaly , There is little day-to-day variation in 
electron densities measured during the summer month's . The electron densities 
measured during the winter months display a great deal more day-to-day 
variability. Enhancements of the electron density above 80 km cause abnormally 
large absorptions of high and medium frequency radio waves on many days during 
the winter months. This phenomenon has been called the winter anomaly. The 
winter anomaly is probably the result of -several causes, and the relative 
importance of each on particular days is unknown at this time. The suggested 
causes of the winter anomaly are summarized by Sechrist [1972J . The causes 
include transport of ionizable minor constituents, the variations in electron 
loss rate associated with positive-ion composition variations, and the effect 
of energetic electrons precipitating from the radiation belts after magnetic 
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storms. Evidence of the winter anomaly in the D- region electron densities 
will be shown in Chapter 5. 

1.2 Electron-Density Measurement Techniques 

The D region is particularly, difficult to obtain measurements from. It 
is too high for balloon measurements, and too low for effective satellite 
measurements. Rocket measurements can be made, but are quite expensive. The 
best alternative for effective long-term study of the region is the use of 
ground-based radio-wave probing. Several radio-wave probing methods are used. 
These include absorption measurements, the ionosonde experiment, the wave- 
interaction experiment, the incoherent-scatter experiment, and the partial- 
reflection experiment. 

The absorption experiment measures the absorption of a radio wave along 
an ionospheric path and can be done in several ways. One method is to trans- 
mit a pulse vertically, which is nearly totally reflected from the E region. 
The absorption due to the D region can thus be determined. The ionosonde 
experiment sweeps through a range of frequencies, transmitting a signal 
vertically. The height where total reflection occurs is measured, and the 
electron density can be calculated. The wave -interact ion experiment [Fejer 3 
1955] uses a high-power transmitter to heat the ionosphere. A second trans- 
mitter is used in probing the heated region. The heating causes a change in 
the electron-neutral collision frequency. The probing signal returns for the 
heated and nonheated ionosphere are then used in calculating the electron 
density. The incoheTent-scatter experiment transmits very high frequency 
signals upwards and measures the scattered energy from the electrons . This 
amplitude can be related to the electron density when a calibration measure- 
ment is used. 

The partial -reflection experiment, which is the subject of this work, 
is described in the next two sections. 

1.2.1 Differential-absorption technique. Ordinary and extraordinary 
mode pulses are alternately transmitted into the ionosphere. These modes 
are the characteristic modes, which by definition propagate through the medium 
without a change in polarization, but undergo a different attenuation. As 
the signals pass through the D region, irregularities cause the signals to 
be partially reflected. By measuring the ratio of the amplitudes of the 
modal returns as a function of height, an electron-density profile can be 
calculated. 


3 , 


1.2.2 Differential-phase technique , In the differential -phase tech- 
nique the phase difference between the modal returns is measured as a function 
of height. The characteristic modes have different phase shifts while pro- 
pagating through the medium. The phase difference as a function of height 
can be used to calculate the electron-density profile. This electron- density 
profile measurement is independent of the differential-absorption measurement. 
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2 . THEORY OF THE PARTIAL-REFLECTION EXPERIMENT 
2 . 1 Generalized Magnetoionic Theory 

The original partial-reflection experiment performed by Gardner and 
Pamsey [1953] used the Appleton-Hartree formula \Budden 3 1961] for the re- 
fractive index of the medium, which assumed that the collision frequency was 
directly proportional to the electron velocity. By laboratory experiments , 
Phelps and Pack [1959] established that the collision frequency of an elec- 
tron in nitrogen gas is proportional to the square of its velocity. Sen and 
WyVlev [1960] derived the formula for the refractive index of the medium using 
a coordinate system with oblique axes, allowing for the harmonic time varia- 
tion of the radio-wave fields with real factors of cos ml and sin ml. Bndden 
[1965] derived the same formula using orthogonal axes and a complex time fac- 
tor giving a simplified approach. The following assumptions are made 

in their derivation: 

1. The medium is a slightly ionized gas. 

2. The mass of neutrals is much greater than the mass of electrons. 

3. The influence of electron-electron collisions is negligible when 
compared with the influence of electron-neutral collisions, and 
collisions are purely elastic. 

4. The collision frequency of electrons is proportional to the square 
of the electron velocity. 

5. The electric field energy is negligible compared with the thermal 
energy so that neutrals and electrons have Maxwell -Boltzmann 
velocity distributions. 

A generalized permittivity tensor is formed in terms of the elements 
£ I* e II* an< ^ e ux» which are the permittivities associated with the principal 
axes defined by the directions of the wave normal and the earth’s magnetic 
field. The permittivity elements are given by: 


I 

= (1 - a) - jb 

(2.1] 

II 

= jCf - d) + (§-](<? - e) 

(2.2] 

III 

= [a - + e)] * 3 '[b - + d ) ] 

(2.3] 
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The semiconductor ‘ integrals 
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M -p- 


v? e~ u du 

2 2 
u - x 


are approximated to within 1% error by the polynomials given by Burke and 
Eara [1963] listed. in Appendix I. 

The generalized refractive index n „ is defined as n = (y - 

O y \C O y uJ 

and is given in terms of the permittivity elements as: 
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where: 


n 


A + B sin^Q + ( B ^ sin^9 - (? cos^9) 


D + E sin 9 


A " 2 e i t®i + e m^ 


B “ e m te i + e m ) * e n 

C = 2GJ e n 
Z> = 2e T 


(2.4) 


E = 2e 


III 


0 = angle between wave normal and the direction of the Earth's 


magnetic field. 

The quasi-longitudinal approximation for small values of 9 greatly simplifies 
the expression: 



where : 


(2.5) 


( 2 . 6 ) 


ta^ = a g cos0 


The quasi- longitudinal approximation is valid for Urbana, Illinois [Pimat 
and Bowhillj 1968] and is used in daily measurements of electron densities. 
The full Sen-Wyller- expression is used in experimental work on the differen- 
tial phase system to avoid the slight error introduced by the quasi- longitu- 
dinal approximation. 

2.2 Differential-Absorption Theory 

The received amplitude of the partially- reflected echo is determined by 
the absorption losses over the path to the reflecting layer and back, and the 
loss on reflection. The two-way absorption is given by [ Budden , 1961] : 
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fh. 

■ 2 L 


(2.7) 


where k is the absorption coefficient defined as ^ x„ „, where y. is the 

O yOO G OyCO 

imaginary part of the refractive index. The received amplitude is then 

/ rh \ 


A « R exp 
o 3 x o 3 x v 


where R is the reflection coefficient. The A JA 
0 3 x' ar c 


( 2 . 8 ) 


ratio can be formed as: 


ar 

iFT exp 


(k x - kj dh 


(2.9) 


Taking the natural logarithm of both sides 


ln r - ln 


¥) * r 

o } } o 


(k - k ) dh 
v x o J 


( 2 . 10 ) 


If the A /A ratio is measured experimentally at two closely spaced 
CO o 

heights h , and h 0> k is approximately constant in that height interval so 

1 b C2jOC 

that the difference in the logarithms of the A /A ratios can be written as 

CO o 


-2(k x - k 0 10t 


C2.ll) 


where Rh - - h^. 

2.2.1 Fresnel reflection model for the differential-absorption 
experiment. The classical partial -reflection theory assumes reflections are - 
due to a sharp discontinuity in the index of refraction. For this Fresnel 
model, the reflection coefficient is given by 


”2 - n l 
n 2 + « x 


( 2 . 12 ) 


Since the observed reflections are weak, the gradient in index of refraction 
is not large and n 2 ~ n l ~ n ~ 1* ^ 0r both ra °des. The reflection coefficient 
can be approximated as: 


o 3 x ~ 2 n 


(2.13) 
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giving us 
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Sn n 
x o 

n 6n 
x o 


(2.14) 


Equation (2.14) is substituted into equation (2.11), with the values for 

the refractive index given earlier, and the experimentally measured values of 

A /A , to obtain electron density values. When using the full Sen-Wyller 
o 

formula for the refractive indices without the quasi-longitudinal approxima- 
tion, there is no closed-form expression for the electron concentration. The 
solution must be obtained by iteration in N. This process requires approxi- 
mately five minutes for the PDP 15/40 computer to compute an electron density 
profile. 

By using the quasi-longitudinal approximation, the electron density can 
be solved for explicitly. 

The values for n are given by equation C2.6), and k is from the 

Oj X 

imaginary part of the refractive index, found by expanding (2.6) by the 
binomial theorem: 
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(2.15) 


(2.16) 


The higher order terms can be neglected. The real and imaginary parts can 

be extracted: The real part, u is 

r *V, x 


o 3 x 


1 - 


The imaginary part, x 



3/2 


'to ± to. 
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(2.17) 


n x 4 oov y? , 
m r~ 


5/2 


'u ± (0, 
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(2.18) 
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S Ne 


in ± to. 


4 m o wv m ^5/2 \ J 

The real part can again be expanded by the binomial theorem, yielding: 


(3.19) 


v e,a “ 1 


0 ) ± 0 ). 


(ii ± 0 ), 


Assuming v is constant across the discontinuity [BeZnose and Burke 3 1964] , 
2 ™ 

n is a function of electron density, B alone. Differentiating equation 
(2.6)' with respect to N yields: 


2n = 
m 6N 


me cov 
o m 


co £ co. 



substituting k „ = ~ with „ from equation (2.19), into equation 

0^30 & 0^30 0^33 

(2.11) results in 


(2.24) 


The collision-frequency profile assumed is given in Appendix II, and 



includes seasonal variations. 


2.2.2 Volume scattering model for the differential- absorption experi- 
ment. The real structure of the irregularities producing partial -reflec- 
tions, is unknown, so there is a great deal of controversy over the type of 
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reflection model to be used. The first approach was the isolated reflector 
[Fresnel) model used by Gardner and Pawsey [1953], and later by Belrose and 
Burke [1964] . Belrose and Burke found that reflections produced by pulses 
of different width showed maximums at the same altitude, with more defined 
minimums displayed with the narrow pulses, and that reflections produced by 
pulses of different frequencies [40 kHz apart) showed similar structures. 
They concluded therefore that the reflections were produced by isolated re- 
flectors. 


Blood [1968, 1969] developed a model with a pulse width dependence. 

His volume scattering model assumes that an echo received at time t is from 
a volume of height ot/ 2, where t is the transmitted pulse width and a is the 
speed of light. He assumes that the return is a result of a volume of dis- 
crete scatterers. The theory allows for the wave to encounter different 
amounts of absorption in the scattering volume. The Flood model can be 
written as [Cohen and Ferraro , 1973] : 



where : 


V 1 - e *P<- 2 X x ■ ,T n {[ « a ,;g 3/2 C««Jl 2 * MW g' 5/2 (<y)] 2 ) 

V 1 - K 3/2 ] 2 * CC5/2) ^ 5/ 2 C“o5] 2} 


• exp 


— 

r°V 2 

-4 

: V* 

4 

o - _ 


( 2 . 25 ) 

* 


K 

K 

XjO 

a 

x,o 


is the average value of K within the scattering volume 
= {[[S/4)ilfe 2 ]/( me Q cv m )} £ 5/2 [0 ± eo^/vj 

= (io ± t 1 i L ')/[v m Cor/2')] 


03? QJ? 

v m (— ) = collision frequency at height — j 
x = pulse width 

r = variable of integration representing time 

= t — t 

t = time measured from the leading edge of transmitted pulse to 

height of interest 

Cohen {1971] and Cohen and Ferraro [1973] give an improved model which 
‘has all the features of the Flood model and. also allows the possibility of 
variation of the reflection coefficient within the scattering volume. The 



expression for the Cohen model is: 
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f r B +r 
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1 2 


ov 
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exp -4 


V s 


where: 



C2.26) 


N (— = electron density at height — 


This reduces to Flood's expression if the electron density and collision fre- 
quency are assumed constant over the pulse width. 

There has been much dispute as to which model is more appropriate \Rolt, 
1969; Floods 1969], and there is no conclusive evidence to support either 
model. Studies of the statistical distributions of the amplitudes have been 
made in an effort to determine the scattering mechanism. It is well known 
that scattering from a single isolated reflector will result in an amplitude 
distribution with a Rician probability density. From the central limit 
theory, scattering from a volume of scatterers would result in an amplitude 
distribution with a Rayleigh probability density. Von Biel [1971] found using 
this analysis that the amplitude distributions were predominantly Rayleigh 
below 80 km, and Rician above 80 km. This would predict scattering due to a 
volume scatterer below 80 km, and scattering from a single reflector above 
80 km. Mathews et oil. [1973] did a slightly more extensive analysis of the 
' amplitude distributions and also predicted volume scattering below 80 km, 
and single reflector scattering above 80 km. Newman [1974] also did an analy- 
sis of amplitude distributions and obtained similar results, predicting the 
scattering mechanism shifting to a single reflector model above 85 km. 

Austin and Mans.on [1969] gave evidence in support, of the single reflec- 
tor model throughout the D region. Further work by Belrose [1970] lead him 
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to conclude that in general the reflection process is not uniquely stratified 
or due to' a volume of turbulence, but is a combination of the two. The alti- 
tudes where each type of reflection occurs varies from day to day, by season, 
and with changes in latitude. 

The differences between electron densities calculated by the volume- 
scattering theory and by the single-reflector theory have been computed by 
Wratt [1974] . He found that for a pulse width of 25 ps and an average elec- 
tron-density profile the difference was at most five percent in the altitude 
range of 70 to 87.5 km. Holt [1969] also investigated the differences in 
electron densities calculated using both reflection models, and found that 
with a 25 ps pulse the difference was less than 10 percent throughout the 
V region. In view of the magnitude of the error that may be introduced, the 
uncertainty of the actual reflection process, and the relative simplicity of 
the single reflector (Fresnel) model, the single reflector model will be used 
in this work,- 

2.3 Differential-Phase Theory 

The received phase of the partially reflected echo is determined by the 
phase change while propagating up to the point of reflection and back, and 
the phase shift due to reflection. The two-way phase shift is given by 
\Budden 3 1961] 

■h • 

p^ dh (2.27) 

o 

where p is the real part of the refractive index of the medium. Defining 

*Xj 

<p.^ and as the phase shifts for each mode due to reflection, we have 



rh 


^o 3 x ~ ^ Ro 3 Rx 


2b) 


p dh 
o 3 x 


C2.28) 


where <J> is the received signal phase. The expression for the phase dif- 

Oj 00 

ference between the two received modes is 


^o' ^Rx ^ Ro 





(2.29) 


Jo 

For closely spaced layers the composition will change only slightly so equa- 
tion (2.29) can be written as a summation over i intervals. 
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2m H a 

k - *o = - ho- — A 

-1 = 0 -lr % 

(2 . 30) 

now defining 



n(&) = - <j> 0 ) 

h 

(2 . 31) 

and writing 



An = n(h 2 ) - n(^) 

(2.32) 

also defining 



= tfite " *Ro 


(2.33) 

and writing 



AB- = B(?j 2 ) - B(^) 

C2.34) 


the contribution to the total phase shift from a slab of thickness A h with 
lower and upper boundaries and respectively is 


An = AB - — ~ (X,, - V 0 )M (2.35) 

2.3.1 Fresnel reflection model for the differential-phase experiment. 
The Fresnel model assumes that the reflections from the ionosphere are due • 
to sharp discontinuities in the index of refraction. The Fresnel reflection 
coefficient can be written as 


R = \R 

OjX 1 c 


R 

e Oj x 


(2.36) 


Thus, the phase, shift due to reflection is the argument of the Fresnel reflec- 
tion coefficient given in equation (2.36). 

Electron-density values are obtained by measuring An, and calculating 

A3 and u in equation (2.35). When using the full Sen-Wyller formula with- 
O 

out the quasi- longitudinal approximation, N cannot be solved for explicitly, 
and the solution must be obtained by iteration in 17. 

The quasi- longitudinaXapproximation greatly simplifies the expressions 
so that the explicit solution for electron density can be obtained. 

From equation (2.21) the reflection coefficients are 
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with the plus sign for the ordinary mode. The argument of R is the phase 
shift Rx and is given by 


*Ro a Rx 

tan -1 



(2.38) 


with the plus sign for the ordinary mode. The expressions for are 

dependent on collision frequency and not on electron concentration. 

The real part of the refractive index, y is given by equation (2.20), 
and from this we obtain 
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0) + u. 


' L m r 3/2 \ m / m r 3/2 V m /J 

Substituting into equation (2.35) gives us an expression for the electron 
density. 
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A priori knowledge of the collision frequency must be assumed. The profile 


used is given in Appendix II. 

2.3.2 Volume scattering model for the differential-phase experiment. 
As discussed in section 2.2.2, the actual mechanism of reflection in the D 


region is unknown. The single- reflector theory for the directly .measured 
differential-phase experiment was given by Connolly [1971] . A volume theory 
following the same basic assumptions as the volume differential -absorption 
experiment of Cohen [1971] was derived by Neuman [1974] . He derived an ex- 
pression for 2?[A<j> ], the expected value of the phase angle contributed by the 
reflection coefficients : 
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where : 


#(-?£) = electron density at 


or 


r or- 


v^C— = collision frequency at — 
= C0} ± “L 3/[v m Cac/2)] 






S 


= variable of integration representing time 
— t - T 

= time measured from the leading edge of the transmitted pulse 
to the height of interest 
= complex random variable geometric scale factor 
= average number of irregularities within the scattering volume 
= &N/N a constant dependence of the irregularities on the back- 
ground electron density. 

The total expected phase difference is then given by 
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Cv x - V 0 )dz + 


C2.42) 


where h is the center of the scattering volume. Newman did an analysis 
of the differences in electron-density calculations for the single- scatterer 
and volume -scatterer models. He found for a 25 ys pulse width that the 
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difference was less than one percent in the altitude region above 70 km. For 
this reason, and due to uncertainties in the actual reflection mechanism, the 
more easily applied single reflector model will be used in this work for 
differential -phase electron-density calculations as well. 
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3. SYSTEM DESCRIPTION 

3.1 System Requirements 

The partial-reflection system is located at the Aeronomy Laboratory 
Field Station near Urbana, Illinois (geographic coordinates: 40° 10' 10" N, 

88° 09' 36" W) . A diagram of the system is shown in Figure 3.1. 

The nature of the signals reflected from £- region irregularities put 
strong requirements on the elements of the system. The reflection coeffi- 
cients are quite small in the D region at the frequency in use, and are of 
the order of 10 -3 at the upper D region and 10~6 i n the lower D region 
[da Silva and Bowhill } 1974] . While the ordinary mode pulse propagates 

i 

through the region, the extraordinary mode pulse is strongly attenuated as 
it passes through ionized layers. At 60 km the reflections from both modes 
are just above the noise level. The ordinary mode returns steadily grow 
stronger with altitude of reflection until they saturate the receiver 
(saturation with full system sensitivity normally occurs at around 85 km) , 
due to the reflection coefficient increasing with altitude. The extraordinary 
mode returns increase with reflection altitude until they peak at around 75 
km, and then at higher altitudes the attenuation effect of the extraordinary 
mode becomes dominant causing a decrease in the amplitude of the returns with 
reflection altitude until they again reach the noise level near 90 km. Also, 
the fading of returns from each height causes the signals to undergo even 
more amplitude variation. These characteristics of the returns set the re- 
quirements for system gain and linearity. 

The system can be broken down into four sections; the transmitter, the 
receiver, the antenna, and the data-acquisition section. 

3.2 Transmitting Equipment 

A large transmitted signal power is required to provide usable returns 
at lower altitudes. The characteristics of the transmitter in use are: 

Peak power: 50 kW 

Frequency; 2,66 MHz 

Pulse width: 23 ps 

Output impedance: 50 OHM, unbalanced 

Pulse repetition rate: 2.5 double-pulse frames/s or 2.5 four-pulse 

frames/s . 

The transmitter is a multistage tube type, and is fully described by Henry 
[1966] and by Pimat and Bowhill [1968] . 
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Figure 3.1. 


Block diagram of the partial -reflection system 
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3.3 Receiving Equipment 

Due to the fact that signal amplitudes must be accurately measured over 
a 50 dB return signal range, an extremely linear receiver must be used. The 
original partial-reflection receiver [Henry 3 1966] had a total range of 20 dB 
for 1 dB deviation in linearity. The detector circuit was modified later 
[Henry in Eduards , 1973a] to achieve a total range of 55 dB for 1 dB deviation 
in linearity. The characteristics of the receiver are listed below: 

Center frequency: 2.66 MHz 

Noise figure: 3 dB 

Bandwidth: 40 kHz (between -3 dB points) 

Recovery time: 200 ps after removal of 0.1 V RMS input 

Gain variation: 3 dB max, 15° to 35° C 

RF input impedance: 50 OHM, unbalanced 

Output impedance: 10 K OHM, unbalanced 

Output response: dc to 50 kHz, 10 V max 

Linearity: 55 dB for 1 dB deviation 

3.3.1 "Phase detection system. This section is a description of the 
phase detector hardware designed by Abert • ' [1975] . To establish 
the proper phase reference, the transmit frequency is mixed down in the same 
manner as the return wave. Figure 3.2 is a block diagram showing how the 
phase reference signal at 455 kHz is derived from the partial -reflection 
system. Using this derived reference signal assures that both received 
return signals and reference signal are equally affected by any possible 
oscillator drift, so as to cancel when their difference is taken in the 
phase detector. 

The mixers used to derive the reference signal are Hewlett-Packard 1053 A 
ring-diode balanced mexers, with a dynamic range of 105 dB. Low-pass filters 
are added after each mixer to attenuate the unwanted mixer products by 40 dB. 
Isolation amplifiers are inserted between the RF and IF oscillators and their 
respective mixer inputs to minimize interaction. The RF local oscillator 
isolation amplifier has an imput impedance of 3 K OHMS, output impedance 
of 50 OHMS and a voltage gain of 2.5. The IF local oscillator isolation 
amplifier has an input impedance of 10 K OHMS, an output impedance of 50 
OHMS, and a voltage gain of 1.0. Figure 3.3 shows the circuitry of the 
isolation amplifiers, mixers, and filters used in deriving the reference 
signal. 
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Figure 3.2 Phase detection system 
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Figure 3.3 Buffer amplifiers and mixing circuitry. 



22 


3.3.2 "Phase detection circuitry t The phase detector must be capable 
d£ working over the large range of amplitudes of returned signals. It must 
also be capable of measuring the range of phase .differences likely to be 
encountered. The expected phase difference has been calculated from an 
electron- density profile obtained by rocket data [Abert in Eduards, 1974b] . 
This ejected phase difference under active sun conditions reaches 180 degrees 
near 83 km and exceeds 360 degrees at nearly 85 km. Under quiet'sun. condi- 
tions, the expected phase difference is below 360 degrees up to 90 km. The 
detector then should have a range of 360 degrees. A phase detector used pre- 

4 

viously by another group [Mathews and ConnaZty , 1972] used a NAND gate as a 
phase comparator, which has a useful range of 180 degrees, leading to ambig- 
uous results at higher altitudes. 

The phase detector in use is shown in the block diagram of Figure 3.4. 

The 455 kHz signal from the receiver IF and the 455 kHz reference signal 
first go through a hard limiter and then a zero-crossing detector to produce 
a square wave output that is independent of input amplitude over a range of 
45 dB. Figure 3.5 is a schematic' of the squaring amplifier and zero-crossing ' 
detector for one of the two channels. 

The phase comparator is a Motorola MC 4044 integrated circuit comparator 
which contains sequential logic circuitry capable of measuring phase dif- 
ferences over a range. of _+ 360 degrees. The circuit of the phase detector 
is shown in. Figure 3.6. The circuit responds only to transitions . and is not 
sensitive to waveform amplitude or 'duty cycle . The output of the comparator 
varies linearly from 0.75 to 2.25 volts as the phase difference varies from 
-360 degrees to +360 degrees. An amplifier stage with level shifting was 
added to make the detector output compatible with the existing analog-to- 
digital converter. This stage provides an output voltage range of 0 to 1 
volt for the phase difference range of -360 degrees to +360 degrees. The 
output amplifier frequency response is flat out to 100 kHz, which is the 
sampling rate of the analog-to- digital converter. A bandstop filter is added 
to the output to eliminate any frequency component at 455 kHz, 

3 . 4 Antenna System 

The partial-reflection antenna system consists of two separate square 
dipole arrays, each made up of 60 half-wave dipoles. One array is used for 
transmitting and the- other is used for receiving, to eliminate the need for 
a transmit-receive switch. Figure 3.7 is a diagram of the arrays. The cal- 
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Figure 3.4 Block diagram of the phase detector 
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culated gain of each array is 22 dB, and the calculated bearawidth is 14 degrees. 
The north-south and east-west set of dipoles are matched to separate feeders 
and brought into the Field Station building. The polarization, for either 
array can be switched inside the building to left-hand or right-hand circular 
polarization, or north-south or east-west linear polarization. Attenuators 
and phase shifters are inserted in each feeder so that the ellipticity of 
each array can be adjusted for the maximum rejection of the opposite circular 
polarization, as measured from ionospheric returns. 

The procedure for adjusting the antenna polarizations -is given in Edwards 
[1973a] . The computer program POLCHK was written to utilize the digital input/ 
output interface to the PDP 15/40 computer to ease the task of measuring the 
rejection of the opposite mode of circular polarization. The program switches 
a manual attenuator into the receiver. antenna input line during the ordinary 
mode return. This attenuator can be switched manually in 1 dB steps to bring 
the ordinary mode return amplitude down to the extraordinary mode return ampli- 
tude, as displayed on the A-scope. Measurements are taken from a sharp dis- 
continuity in the E or F region. Thus, the rejection ratio of each' antenna 
system can be determined. The program includes the provision for setting the 
digital attenuator of the system to any desired attenuation up to 63 dB, to 
prevent saturation of the receiver from these strong signals. To change the 
attenuator setting, data switch 17 depressed momentarily and the desired at- 
tenuation is typed into the teletype followed by a carraige return, - The pro- 
gram has been added to the main partial-reflection program tape, and all of 
its subroutines (INTER, INPADF, SYNC, OUT) are included in the library file 
on that tape. To run the program, .DAT slots -4 and -5 are assigned to the 
tape, and the linking loader is called to load and execute POLCHK. A list- 
ing of the program and its subroutines is given in Appendix III. 

3.5 Data^Acquisition System 

A Digital Equipment Corporation PDP 15/40 digital computer is used for 
data acquisition and processing. It is a background-foreground system with 
18-bit words and 32 k of core memory. Peripherals include two disk pairs, 
four DECtape drives, a paper punch and reader, a cathode-ray- tube display 
terminal, and a General Electric Terminet-1200 high-speed teletype. A 
Hewlett-Packard 5610 A analog- to- digital converter is used to digitize data. 

It has a multiplexer so that any one of 16 channels can be sampled at a time, 
and has a 10-bit word output, in two’s complement form. Thus, the 0 to 1 
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volt input range in use corresponds to an A/D output of from 0 to 512 base 10. 
Conversion of data in the free run mode occurs at a 100 kHz rate, corresponding 
to an apparent height resolution of 1.5 km, for a single channel input. A 
complete .description of the computer system can be found by Birley and Seehris t 
[1971], and Bean and Bowhill [1973]. 

3.5.1 Sample-and-hold circuitry for collection of simultaneous amplitude 
and phase data ; A circuit has been built to allow simultaneous amplitude and 
phase data collection. The circuit shown in Figure 3.8 contains two sample- 
and-ho'ld circuits, and the associated timing circuitry. The Hewlett-Packard 
5610 A analog-to-digital converter has a conversion time of 10 ys. It is used 
in the free run mode, so a sample is converted every. 1.5 km of apparent height. 
For simultaneous operation, the external sample-and-hold circuits will be trig- 
gered every 20 ys (or 3 km of apparent height) , in the desired height region. 
These' two inputs will be held on the A/D input channels 0 and 1. With the 
A/D converter in the free run mode, an encode command, given in the height 
region 45 to 90 km, causes samples to be converted every 10 ys. The channel 
sequence logic, itfhich externally selects which of the 16 channel inputs on 
the A/D is to be sampled [Hess in Edwards, 1975a; Hess and Bcwhill, 1976] was 
modified to correctly sequence through the two input channels. Figure 3.9 
shows the inputs added to the channel sequence logic. The reset to channel 
0 assures that collection will start with- amplitude data at the first mea- 
surement. The. data collection program will sort out the sequentially con- 
verted samples for the processing programs. 

\ 

The sample-and-hold command and channel sequencing signals are derived 
from the data-ready flag of the A/D converter. Therefore, the circuit is 
active only during the encode command period. A timing diagram is shown in 
Figure 3.10. The chain of monostable multivibrators is used to obtain a 
1-ys-wide pulse delayed from the data-ready flag by 4.5 ys. A J-K flip-flop 
is used to divide the number of delayed pulses to provide a single pulse 
every 20 ys. The flip-flop is preset to the 1 state by a transmit time 
pulse, so that the first data-ready flag will be followed by a sample com- 
mand during the next conversion cycle, and conversion of the inputs will 
occur during the second next two conversion cycles, To obtain the correct 
height for the samples, the .encode signal must be given to the A/D converter 
14 ys earlier than in the present system (which starts conversion 2 ys after 
the encode command) , so that the first sample will occur at the desired ap- 
parent height of 45 km. 
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Figure '3.9 Modified channel sequence logic 
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4. DATA- ANALYSIS TECHNIQUES 
4.1 Analysis of Differential-Absorption Data 

Differential- absorption data are collected on a regular basis at the 

Aeronomy Laboratory Field Station. Data are collected daily near noon during 

/ ' 

the months of October through April, and twice weekly during the summer months 
These data are being used in studies of day-to-day D - region electron- density 
variability associated with the 'winter anomaly'. For these studies, a one- 
hour collection routine is used to obtain a value of electron density at 72, 
76.5, and 81 km, and a value of A X /A Q at 81 km. These data are plotted and 
published semiannually [Edwards^ 1973a, b, 1974a,b, 1975a,b, 1976a, b, 1977a,b, 
1978] . The collection routine summarized here is described in detail by Bean 
and Bowhill [1973] . * 

The transmitted frame consists of a pair of pulses, an ordinary mode 
pulse, followed 33 ms later by an extraordinary mode pulse. Data are digiti- 
zed by 1.5 km intervals over the height range of 45 to 90 km, and collected 

in 3.4-minute' files . The data from the lowest five altitudes, from 45 to 

> 

51 km (where no appreciable electron density exists, and hence no returns 
are expected] are used to estimate the noise level during the collection of 
that set of data. A noise threshold is determined from the previously col- 
lected 3.4-minute file of noise samples (the first file uses noise data from 

45 samples taken initially) . The square of this noise threshold is multiplied 

, < 

by five and compared with the sum of the squares of the five low-altitude 
noise samples from each data set. If the noise of either the ordinary or the 
extraordinary mode is greater than the noise threshold, both sets of data 
are rejected, A running count is kept of the number of frames rejected, to 
provide an indication of the actual number of frames of data used in each ' 
file, The noise power obtained from the noise samples is then subtracted 
from the data returns. Data digitized in the height range of 52.5 through 
58.5 km are not retained for further processing. 

Files of 513 samples (or 3.4 minutes) are collected with 0 dB, 10 dB, 
and 25 dB of attenuation inserted ahead of the receiver. .An electron-density 
profile is calculated for each file. This group of -three files is merged to 
obtain a single electron- density profile, by receiver saturation criterion. 

The receiver saturates at one volt output, so any data above this amplitude 
is assumed to have saturated the receiver. A running count is kept of the 
number of saturations at each altitude. If more than ten samples saturated 
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the receiver at a particular altitude, the data at that altitude are then 
taken from the next higher attenuation file to form the merged file. 

Five of these merged files are collected, requiring about one hour. The 
median value of electron density is found at each altitude from the five 
merged files, and a final electron-density profile is printed. Averages are 
calculated for 3-km-wide slabs centered at 72, 76.5, and 81 km. This is 
done by averaging the electron- density values from two adjacent altitudes. 

For example, the 72-km electron density is obtained by averaging the values 
from 70.5 and 72 km. Since the electron density at 70.5 is obtained from 


ln( A x /A q ) 


~ ln( VV 

72 km 


70,5 km 


(4.1) 


and the 72-km electron density is obtained from 


lft (Ac/ A?) 


- In (A^./ A^) 

73,5 km 


72 km 


(4.2) 


The average of these two is representative for 72 km. These average values 

* 

at 72, 76.5, and 81 km are plotted for further study, as is done in Chapter 
5. 

4 . 2 Analysis of Differential-Phase Data 

This section deals with the characteristics of the phase detector returns, 
and the methods used in collection and processing of the data. 

4.2,1 Characteristics of the phase-detector returns . The output of the 
phase detector is a phase angle with respect to the system reference oscil- 
lator, The phase angles measured have an ambiguity of 2m (n an integer), so 
that the angles measured at two adjacent altitudes may be written as 

$ + 2air (4.3) 

*i 

+ 2bir . (4.4) 

°1 

<j> + 2eir (4.5) 

<J> + 2dir (4.6) 

°2 

x, o denoting mode; numerical subscript denoting altitude; a 3 b 3 c, d are 
arbitrary integers . 

Consider from equations (2.31) and (2.32) 

An = (^2 " " ^0]_) 


( 4 . 7 ) 
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where <j> , <j> > , and are the true values of the phases . From cal- 

x 2 

culations using electron-density profiles obtained from rocket data [Edwards 3 
1974b (Table 6.1)-] it is- expected that An will- lie between -ir and +ir up to 
about 83 km under active sun conditions and 85 km under quiet sun conditions, 

i 

with h 2 - hi - K5 km. Thus, a multiple of 2ir is added to the phase dif- 
ference An to obtain a value between -ir and +ir. 

Figure 4.1 shows two typical histograms - of phase difference An,’ at 73.5 
to 75 km, and 78 to 79.5 km. The data were collected! near noon with a sample 
size of 5000'. (The systematic variation between adjacent bins is due to the 
quantization error in the A/D converter and histogram bins.) 

Histograms obtained from altitudes above 80 km are much broader than 
those displayed in Figure 4.1. This leads to an ambiguous determination of 
the phase angle, leading to large errors in the electron- density calculations 
at those altitudes. Similar results were reported by Austin [1971], He at- 
tributed this effect to the different group velocities of the two propagating 
modes. The group velocities increasingly diverge as the difference between 
their refractive indices becomes larger with altitude. Austin defines an 
overlap function to describe the degree of overlap of the scattering volumes 
of the two modes. He finds that due to the group delay there is total over- 
lap of scattering volumes- up to about 86 km, and no overlap above 9$ km. 
Tanenbaum et at. [1973a, b] made a study of phase -difference histograms and 
tried to fit them to histograms predicted by single-reflector and turbulence- 
volume models. They found that the observed data were in good agreement with 
both models,, and could not predict the type of scattering. They also obtain- 
ed phase-difference histograms which broadened with increasing altitude. 

They concluded that the broadening of the higher altitude phase-difference 
histograms was not necessarily due to nonoverlap of the scattering volumes, 
but due simply to the increasing difference between the differential phases 
of signals scattered from the top and bottom of the scattering volume. 

Due to this'broadening of the phase-difference histograms, reliable 
electron densities cannot be obtained above about 80 km. Two methods of 
extending the altitude range of the experiment would be to use a shorter 

pulse length or a higher operating frequency, both of which would limit the 

\ 

phase angle difference to smaller values. 

4.2.2 Data collection and 'processing programs. The computer program 
used in collection of the differential -phase data is PCOLL. Listings of the 
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program and all its subroutines are given in Appendix III. The subroutines 
are included in the partial-reflection program library C.LIBR5) and may be ^ 
loaded automatically by the linking loader. 

The program collects the data from the A/D converter and writes it in 

files onto a magnetic disk or DECtape. Provision is included to switch the 

system digital attenuator in to put any attenuation up to 63 dB ahead of the 

* 

receiver to prevent its saturation. Normally, differential-absorption data 
are obtained first to obtain the attenuation necessary to prevent receiver 
saturation in the altitude range of .interest . Data are stored in files of 
any length, as selected by the operator. 

The phase angle processing and plotting program is PHIST, and is listed 
in Appendix III. The program reads the file of data written by PCOLL, cal- 
culates the phase difference An, adds or subtracts 2ir to bring the phase 
angle into the range of -ir to +ir, and plots the results in 1.5-degree incre- 
ments . 

The differential -phase electron- density profile calculation program is 
DPEX. This program and its many subroutines are also listed in Appendix III. 
DPEX reads the data file written by PCOLL and calculates the electron-density 
profile in the altitude range of 70.5 to 90 km. The program obtains a histo- 
gram of the phase difference An at each pair of heights and finds the median 
of the distribution. This is called the first median. A 'second median' 
is obtained by taking the median of the distribution truncated at plus or 
minus ,20 degrees from the first median. This is the phase difference An used 
in the electron-density profile calculation. This 'second median' is used 
because it was found through experiment to be the best approximation to the 
center of the distribution of those methods tried. The electron densities 
are calculated using the full Sen-Wyller theory, by iteration in electron 
density. The profile calculation requires about five minutes of computation 
on the PDP-15/40 computer. 
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5 . PARTIAL-REFLECTION ELECTRON-DENSITY MEASUREMENTS 

In this- chapter electron-density data analyzed by the methods described 
in Chapter 4 are presented . 

5.1 Daily Measurements of Electron Densities 

In order to identify the changes that occur in the Z?-region electron 
densities, data are collected on a daily basis at the Aeronomy Laboratory 
Field Station during the months of October through mid-April. Figures 5.1 
through 5.7 are plots of the electron density at 72, 76.5, and 81 km, and 
the A X /A Q ratio at 81 km for the months of October 1975 through April 1976. 

The A x /A 0 ratio, whose logarithm represents the total differential-absorption 
occurring below 81 km, gives an indication of the total electron density below 
81 km. A high value of A x /A 0 at- 81 km corresponds to a low total electron 
content below 81 km. The variations in electron density are relatively minor 
until the second week in November. Variability is greatly increased there, 
and the variability continues to grow larger until it maximizes in the month 
of February. This is clear evidence of the winter anomaly occurring at 
Urbana. Variability is again lowered by the second week in March, and con- 
tinues at a low level through mid-April . When comparing plots of a particular 
altitude with plots of another altitude it can be seen that the electron 
densities at 72 and 76.5 km follow the same variational trends, suggesting a 
similar mechanism causing the changes at both altitudes. The 81-km electron 
densities are not as well correlated with those at the lower altitudes, sug- 
gesting that different mechanisms are causing the variation at 81 km. It 
can be seen though that when there is a sizeable peaking of electron density, 
it usually occurs at all altitudes. This suggests a mechanism that is large 
enough in consequence to disrupt the normal variations observed at each al- 
titude . 

Figure 5.8 is a plot of the electron density monthly means at 76.5 km. 

The electron densities are shown to be higher in the summer months as would 
be expected due to a smaller solar zenith angle in that period. Error bars 
indicating +_ 1 standard deviation demonstrate the increased winter variability. 
The coefficient of variation (standard deviation/mean) is plotted in Figure 
5.9 for electron densities at 72, 76.5, and 81 km. The increase in vari- 
ability during the winter months associated with the winter anomaly is 
clearly shown. Gregory and Manson [1969b] did a similar analysis in the 
Southern Hemisphere and obtained results in agreement with those presented 
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lot of electron density at 72, 76.5,- and 
id A x /A 0 at 81 km for December 1975. 
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Figure 5.4 Daily plot of electron density at 72, 76.5, and 
81 km and A x /A 0 at 81 km fox January 1976. 
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FEBRUARY 

Figure 5.5 Daily plot of electron density at 72, 76. S, and 
81 km and Ax/A 0 at 81 km for February 1976. 
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Figure 5.7 Dg f electron density at 72, 76.5, and 81 km 

and A x /A q at 81 km for April 1976. 
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Figure 5.8 76.5-km electron-density monthly means. Error bars 

indicate + 1 standard deviation. 
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Figure 5,9 Coefficient of variation of electron densities at 72.0, 76.5, 
and 81.0 km. 
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